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QED3 with fermi and gap anisotropies is considered to be a candidate effective field theory of 
high temperature superconductors. Simulations of a variant of the theory have demonstrated that 
there is evidence consistent with a phase transition to a chirally restored phase as the velocity 
anisotropy increases, and that the correlation lengths of pions in spatial directions behave in a 
manner consistent with the anisotropy being a relevant parameter of the theory. We present the 
first measurements of the fermion propagator in Landau gauge for this theory and suggest that the 
structure of this theory is more complex than previously thought. 
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1. Introduction 

QED3 has been the subject of recent interest in condensed matter theory as a candidate effec- 
tive theory of the underdoped phases of cuprate high-temperature superconductors [|l], 0J. 

These effective theories have an interesting peculiarity: since we have no reason to suspect 
that Lorentz invariance is respected in a condensed matter system, we can introduce 'velocities' vf 
and va that parametrise the deformation of the Fermi surface of the underlying condensed matter 
theory. At low energies, our action is [Qj: 

f 3 
S= J d 2 rdT{ x i , i[ro(d t + ia t ) + -y=Yi(d x + ia x ) + 8y/Kj 2 {d y + ia y )] x ¥ 1 



+ {\^2,x^y) + -^F llv F llv } (1.1) 



where is an excitation around one of the two pairs of nodes possessed by the d-wave Fermi 
surface, K = ^ (this quantity is proportional to the doping of a cuprate sample), 8 = y/v^Vf and 
the gauge fields model statistical interactions between excitations at the nodes. 

In this approach, it has been suggested that the chirally-broken phase of the model corresponds 
to the onset of a 'spin-density wave' phase ^ [|], and that the chirally symmetric phase to the so- 
called 'pseudogap' phase, some of whose peculiar properties, it is hoped, may be explained by 
chirally-symmetric QED3 [jlj|. 

In order that we might make use of an exact HMC algorithm, we use the following simplified 
version of dO| ), which we expect to capture most of the physics of the system: 

where M x x i is the staggered fermion matrix, modified as follows, 

I 3 

Mxj = — £ ZpWiSjt f+pJJxn - d^^U], ) +m8 tiv , (1.3) 



where ^(x) = X^r\^{x): tj^(x) = (— l^i+'-'+^-i, which is the Kawamoto-Smit phase factor ( 
where x\ =x,X2=y and X3 = t), and A* ; = , Xy = K~- , Aj = 1. 

U X fi = exp(ia6 x u) is the gauge parallel transporter, and the gauge action has the non-compact 

form 

MV (x) = i[A+0 v (x)-A+0 /i (x)]. (1.4) 

In we presented evidence from simulations on a 16 3 lattice of a phase transition to what appears 
to be a chirally symmetric state as K increases above a critical value, fQ. « 4.35, and that, contrary to 
previous analytical studies, the renormalised value of K (when defined as the ratio between the pion 
screening lengths in the x- and y-directions) is a relevant parameter (i.e. (K nren — 1) > (»Q, are — 1))- 
In what follows, we present measurements of the fermion propagator in Landau gauge, whose 
behaviour is thought (not uncontroversially) [jl|, [5j], to provide some explanation for the behaviour 
of the pseudogap phase. 
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2. Gauge Fixing 

The fermion propagator is not a gauge invariant quantity; prior to measuring it we must nu- 
merically impose a gauge fixing condition. We do so as follows [f7|, ||] : 

• We impose 'minimal' Landau gauge, corresponding to the continuum condition d^d^ = 0: 

This corresponds to the minimum of a functional F[6]. For reasons of practicality, in numer- 
ical simulations one imposes the condition up to a residual value, R, which we choose to be 
the floating point value 10~ 6 . 

• To remove all remaining gauge degrees of freedom (largely due to zero modes that are an 
artefact of the finite volume of the lattice), we rotate 6^ x i— ► 6^ x ±nj^, VQ^^/i, where n is 



an integer and is the size of the lattice in that direction, such that: 



2n - 2n 

-r<Qix< t, (2-2) 



where = \ Zx^x ■ 



We call this the 'modified non-compact iterative Landau gauge' or 'miLandau gauge' (by 
analogy with the similarly named condition proposed in [^] for compact QED). 

It is, of course, necessary to check that there are no false minima of F[6] (Gribov copies) that 
might contaminate our measurements; while this gauge condition has been used before, this check 
has yet to be carried out for noncompact QED3, simply as it seems unlikely that there could be any 
such problem. It would be good to confirm such an intuition (particularly as compact QED3 has a 



non-trivial Gribov problem QlOp ) 



In order to do this, we followed the method of [11]: From a set of mother configurations, we 



produced 3 sets of 500 daughter configurations via the following gauge transformations: 
A: Bfix l— ^ @u.x ~ Ohc+fi + OCx, where a x is a random number. 

B: 0fx x 1 — t d^ x ± nj^-, MB^x in the direction ji, where n is randomly chosen to be 0, 1, or -1. 
C: We perform transformation A and then B. 

We then define a 'variance' dF = max !; - [F [{?],- — where i, j run over a mother configu- 

ration and all its daughters. If, for a given mother and its offspring, dF < R, then for all practical 
purposes, there are no Gribov copies. Our results are presented in Figure |l| for a 16 3 lattice with 
m = 0.03, j8 = 0.2 and K = 1.00 (we obtained similar results for a value of K = 10.00). We see that 
for set A, there are no gauge copies for either gauge condition. However, for B and C, we see that 
there are Gribov copies present for mLandau gauge, but not for miLandau gauge. Furthermore, the 
variance is more-or-less identical for both of these groups, which implies that the source of Gribov 
copies are rotations of kind found in group B. We conclude that, in line with our expectations, the 
miLandau conditions fix the gauge completely. 
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Figure 1: Variances for the sets of gauge transformations A, B and C. 

3. Fermion Propagator 

Fermion propagators for j8 = 0.2, m = 0.03 on a 16 3 lattice were measured over 30,000 tra- 
jectories for 1 < K < 10, the extremely large number of trajectories being necessitated by the 
noisiness of the observable. A very slight tendency towards sawtooth behaviour was noted in the 
time-direction propagators above K c as the value of the propagator on even-numbered slices began 
to drop; this was more pronounced in the y-direction and absent in the A-direction. It is not clear 
why this should be; it may have something to do with the restoration of chiral symmetry, though it 
does not quite replicate the expected behaviour of a fermion propagator under such circumstances 
(JT2|] gives an example of what we expect) . Ideally, one would like to fit data with this behaviour 
to a four parameter fit; however, these fits proved to be unstable on this lattice size, and so we fitted 
a two-parameter curve to the odd timeslices only. 

Our results are presented in Figures || and |3[ where fCyy, the renormalised anisotropy with 
respect to fermions, is the ratio of the fermion screening mass in the A-direction with that in the y- 
direction. We may draw two conclusions here: firstly, that fermions retain a dynamically generated 
mass above fc c , and secondly that Kf ren — 1 < Kb are — 1 - that is, K is an irrelevant parameter with 
respect to fermions. This last is quite puzzling, given that K Mren is a relevant parameter of the theory, 
and may take some thought to resolve. 

This retention of the dynamically generated mass into an apparently chirally-symmetric phase 
is more readily explicable, however, and has also been observed in the 2 + ID Gross-Neveu model 
at finite temperature ]T2|]. Witten [13] suggests that such behaviour may be observed if: 



The physical fermion corresponds to a branch-cut, and not a pole, and is constantly interact- 



249/4 



Fermion propagators in QEDt, with velocity anisotropics 



Iorwerth Owain Thomas 



1.5 



0.5 



[♦ mass=0.03 



2 4 6 8 10 

K 

Figure 2: Fermion masses for various values of K. 
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Figure 3: The renormalised anisotropy, Kf rm , for various K. 
ing with massless, non-Goldstone bosons. 

(YY}' tne chiral condensate, is close to having long-range order - for example, it has a phase 
parameter that slowly varies with spatial/temporal position. This would entail an overall non- 
zero magnitude of the order parameter, but fermions would still experience a dynamically 
generated mass. 



If correct, this behaviour is unanticipated by any analytic prediction that we are aware of. 
Because of this, it renders the status of the condensed matter model a little unclear. 

Nevertheless, more research is needed before this can be clarified - the thermodynamic and 
continuum limits must be taken, for example, and it would be good to have some measurements of 
the fermion mass closer to the chiral limit. Even if this theory proves not to be viable as a model of 
a condensed matter theory, we hope that its behaviour is interesting enough in its own right to be 
deserving of study. 
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